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THEORY OF THE TWO-SUBLAHICE HEISENBERG ANTIFERROMAGNET 
by Edwin G. Wintucky and Larry Flax 
Lewis Research Center 

SUMMARY 

Thermodynamic properties of the two-sublattice anisotropic Heisenberg antiferro- 
magnet are calculated by the Green function method in the random -phase (Tyablikov) and 
Callen approximations for spin 1/2. Analytic expressions for the sublattice magneti- 
zation and internal energy are obtained which are valid over the whole temperature range 
up to the Neel temperature for the body -centered cubic lattice. 


INTRODUCTiON 

Antiferromagnetism, like ferromagnetism, is a cooperative phenomenon with an 
ordered state below a well-defined critical temperature. For the antiferromagnet, how- 
ever, the exchange energy is negative, and neighboring spins tend to aline antiparallel. 

The simplest and most often used model for the antiferromagnet is the two -sublattice 
model wherein the lattice of magnetic atoms is subdivided into two equivalent interpene- 
trating sublattices. If only nearest-neighbor interactions are considered, the spins on 
one sublattice interact only with the nearest spins on the adjacent sublattice. A negative 
exchange interaction energy results in the two sublattices being spontaneously magnetized 
in opposite directions . This is the model originally proposed by Neel and has been much 
used in subsequent theoretical studies of the antiferromagnet. At absolute zero the spins 
on each sublattice are completely alined antiparallel, and with increasing temperature 
the spontaneous magnetitization of each sublattice decreases and finally disappears, in 
zero applied field, at a critical temperature called the Neel temperature. Because of 
the equivalence of the two sublattices, the net magnetization is zero. The two -sublattice 
model, which is considered in this report, is compatible only with simple cubic and 
body -centered cubic crystal structures. More complicated models must be used for 
antiferromagnetic materials with the face-centered cubic or other crystal structures. 
Ftu-ther discussion of the two-sublattice and other models, as well as the experimental 


properties of antiferromagnets , can be foimd in the books on magnetism by Morrish 
(ref. 1) and Smart (ref. 2). 

The theory of antiferromagnetism has been much investigated in recent years. The 
Green function technique is now well established as one of the most useful statistical 
methods in the theory of magnetism and has been applied many times already to the study 
of antiferromagnetism (ref. 3 and references cited therein). Earlier methods, espe- 
cially the molecular field theory, are discussed in detail in references 1 and 2. A use- 
ful featiire of the Green function method is its applicability over the whole temperature 
range. Calculations of thermodynamic properties using Green functions have been found 
to be in good agreement with the results of noninteracting spin wave approximations at 
low temperatures and with the results of other statistical methods near the Neel tem- 
perature (ref. 3). 

In this report the Green function technique is applied to the anisotropic Heisenberg 
exchange model of a two -sublattice antiferromagnet to obtain the correlations between 
spins, from which the thermodynamic properties can be calculated. 

The Heisenberg model is considered reasonably valid for insulators and semicon- 
ductors. The sublattice magnetization and other thermodynamic quantities involve com- 
plicated lattice sums which, generally, have been evaluated by means of series expan- 
sions at the high- and low-temperature limits. In a few cases numerical calculations 
have been made at a limited number of intermediate points (ref. 3). Flax and Raich 
have developed a method for evaluating the lattice sums and successfully applied it to the 
ferromagnet to obtain analytic expressions for the spontaneous magnetization valid over 
the whole temperature range (refs. 4 to 6). In this report, their method is applied to 
the antiferromagnet with spin 1/2 to obtain analytic expressions for the sublattice mag- 
netization and internal energy for the body-centered cubic (bcc) structure. The calcula- 
tions are made for two Green function decoupling schemes, the random-phase approxi- 
mation (RPA), which is the Tyablikov approximation, (ref. 7), and the Callen approxi- 
mation. The random-phase approximation is justifiable at low temperatures, where the 
spins are nearly completely alined and fluctuations are indeed small. However, near the 
critical point, fluctuations in become significant, and the RPA gives an overestimate 
of the transition temperature. The Callen approximation takes these fluctuations into 
account and, for spin greater than 1/2, gives better results near the critical point. The 
Neel temperatiures are also calculated and compared with previous results. 


GREEN FUNCTION METHOD 

As stated previously, the model used in the present study is a cubic, two -sublattice 
anisotropic Heisenberg antiferromagnet with nearest -neighbor interactions. The 
Hamiltonian is 
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where ^ is the exchange energy, s'^ and S* are spin operators, and al and 
refer to lattice sites I and m in sublattices a and /3, respectively. Here 
a, ^ = 1, 2 with N/2 sites in each sublattice', t] is an anisotropy factor where tj = 0 
corresponds to the Ising model and 77 = 1 to the isotropic Heisenberg model. In the 
Green function method, the sublattice magnetization and internal energy <H) are 
calculated from correlation functions involving pairs of spin operators, for example. 

For a general description of the double -time temperature dependent Green 
fimction method, see Zubarev (ref. 8 and references cited therein). (Symbols are de- 
fined in appendix A . ) 

The eqmtion of motion in the energy (E) representation for the Green function 
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The commutators for spin operators 
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where 6^^ and 6^^ are Kronecker deltas, are used to evaluate , hJ . The equa- 
tion of motion is then 
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(4) 


Equation (4) represents an infinite hierarchy of coupled equations since each Green 
function on the right side is of higher order and satisfies a similar equation of motion. 
To solve this system of coupled equations some decoupling approximation must be made 
to limit the number of equations. We treat the two decoupling schemes most commonly 
used in the theory of ferro- and antiferromagnetism, the random-phase and Callen ap- 
proximations . 

Decoupling the equations of motion enables one to solve for the Green function in 
terms of the energy and excitation energy spectrum. The corresponding correlation 
fimctions are related to the Green functions by the following eqmtion: 


<B(t’)A(t)> 
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( 5 ) 


where A and B are the appropriate spin operators, t and t' are times, /3 = lAT, 
k is Boltzmann's constant, and T is absolute temperature. Thermodynamic quantities 
of interest, for example, sublattice magnetization and internal energy, are obtained 
from these correlation functions. 

SUBLATTICE MAGNETIZATION 
Random -Phase Approximation 

In the RPA, the fluctuations in are ignored and the operator is replaced by its 
average value. Then 
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The equation of motion becomes 
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where 


^al , |3m 
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Assumir^ translational invariance of the direct lattice, we Fourier transform to the re- 
ciprocal lattice. Let 


^al 
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where k is a vector in the reciprocal lattice, and R^^ is a vector in the direct lattice 
for lattice site I in sublattice a. Fiurthermore, we make the nearest-neighbor ap- 
proximation in which the spins on sublattice 1 interact only with spins on sublattice 2. 
Let 


A sR -R 

al ng 


( 11 ) 


denote a nearest -neighbor vector. Then 

J 7 = J6 < 6 7 A 
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( 12 ) 
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Therefore , 
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(13) 


where J(0) = Jz, 


y(k) 


;Z-“' 


(14) 


is the structure factor, and z is the number of nearest neighbors. The equation of mo- 
tion becomes 


E + A^^(k) = J(0)y(k)A^^(k) {a = a ± 1) (15) 


or 


E + K 


l(s?)] = i(s"^ 6^^ * K2(s^ r(k)A^^(f) 


(16) 


where = 2J(0) and Kg = 2J(0 )tj. For S = 1/2 the sublattice magnetization can be 


written as m = /Sf\ = -/s»\ . 


Substitution into the equation of motion yields two pairs of simultaneous linear 
equations . One pair gives the Green functions 
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where 


Ej^ = m 
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is the excitation energy spectrum. The Green fimctions Agg and are similarly 

obtained but are not needed in the calculation of the thermodynamic quantities , 

The corresponding correlation functions are 
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= m(-l + 
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Transforming back to spin operators, we get 
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and 
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( 23 ) 



The lattice sums appearing in equations (22) and (23) are taken over the first Brillouin 
zone and can be evaluated analytically by the Flax-Raich method. The results are ex- 
pressions for the correlation functions which are valid over the whole temperature range 
up to the Neel temperature. 

The series expansion 


coth(Trt) = 



R=1 


R2^t2 


(24) 


is used where Trt = ^Ej^/2. For a cubic two-sublattice model 


N 



. . . dx dy dz 


The correlation fimctions can then be written as the sums of integrals: 

= m(-l + Ij + I 2 ) 


(25) 


and 


(s;s^ = -Wj * n^) 


(26) 


where the integrals Ig, and Ilg are given in appendix B by equations (Bl) to (B4). 
When evaluated for the bcc lattice, the integrals become 
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where 



and K(/c) is a complete elliptic integral of the first kind. For spin S = 1/2 the sub- 
lattice magnetization is 


" (Si) 4 ■ 


Then 


m = 


1 
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The RPA sublattice magnetization m is plotted as a function of reduced temperature 
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T = kT/J(0) in figure 1 for 77 = 0, 0. 5, and 1. In the RPA, the correlation fimction 

is not needed in the calculation of the sublattice magnetization but is used later 
for the internal energy. 


Callen Approximation 

The Callen approximation is a symmetric decoupling scheme which attempts to take 
into account the fluctuations of around its average value and was first proposed for 
the ferromagnet (ref. 9) and subsequently applied to the antiferromagnet (refs. 3 and 10). 
For the antiferromagnet, 




For 


In the case of the ferromagnet, the empirical factor a is chosen to be 
the antiferromagnet, two choices for a are possible because of the existence of the two 
sublattices, namely and the absence of an applied field both 

choices lead to identical results. In the present calculation, we will use 
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After decoupling and Fourier transforming to reciprocal space, the equation of motion is 




a0 ■ {^ 1 “) 










J (k-k’)!// (k’) 


iix / a) 


(38) 


10 



Applying the same nearest-neighbor approximation used in the RPA calcxilation gives 


EA^^(k)=i(s^)6^^-2(sJ>J(0) 


NS 
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where (ref. 6) 

^ J(k - = J(0)y(k)^ y(k*)l//^^(k') 

k’ k’ 

By symmetry, 

^^2^k) = 

Let 

= 2J(0)(1 - 2t?F) 

Kg = 2J(0 )(tj - 2F) 

where 

k’ 


and S = 1/2. Then 

(e . K, (s^)) A^p(ir) = i(s^a„p * Kj(s^)v(k)A^^(f)(M ^ a) 


(40) 


(41) 

(42) 

(43) 

(44) 
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The equation of motion for the Callen approximation has the same general form as for 
the RPA and leads to a similar set of coupled equations for the Green fimctions . The 
correlation functions are given by equations (20) and (21) . 

The difference is that for the Callen approximation Kj and Kg are defined by 

equations (42) and (43), respectively. The second correlation fvmction is iden- 

tical to the sum F which appears in and Kg and is given by eqmtion (47) . In the 
same way as for the RPA, the correlation functions in the Callen approximation can be 
written as the sums of integrals as in equations (25) and (26). The integrals have the 
same form for both the RPA and Callen approximations and are given in appendix B by 
equations (Bl) to (B4) . The results are given by equations (27) to (30) with K^^ and Kg 
defined by equations (41) and (42). Also, 


(46) 



where 





(47) 


Q := H (1 - 27?F) (48) 

T 

and T is defined by equation (31) . The sublattice magnetization m for the Callen ap- 
proximation is plotted against t in figure 2 for ?? = 0, 0. 5, and 1. A comparison with 
previous theoretical calculations (ref. 3) shows that for spin 1/2 the RPA gives a better 
approximation of the sub lattice magnetization. 


INTERNAL ENERGY 

The internal energy is a useful thermodynamic fimction for the calculation of many 
other properties. For a system of interacting spins, differentiation with respect to tem- 
perature gives the specific heat due to spin disorder. The entropy and other properties 
may then be calculated. In previous studies of the antiferromagnet, the internal energy 
has been derived only for very low temperatures or for the paramagnetic state above the 
Neel temperature. In this section an exact expression is derived for the internal energy 
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of a two-sublattice anisotropic Heisenberg antiferromagnet with spin S = 1/2. This ex- 
pression for the internal energy is then evaluated in the random-phase and Callen ap- 
proximations for the bcc lattice by means of the same method used for evaluating the 
lattice sirnis in calculating the magnetization. The resvilt is an analytic expression for 
the internal energy valid over the whole temperature range up to the Neel temperature. 
The calculation is made for zero applied magnetic field, but the method is equally useful 
for nonzero fields . 

The internal energy U is obtained from the average value of the interacting spin 
Hamiltonian; 


U=<H> 


' ((®lf®2g) ^ (^1^2g) 


(49) 


The correlation function F = has already been calculated, and 

=NJ(0)F 
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is calculated exactly for spin 1/2 in appendix C and is 
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An exact expression for the internal energy of the two-sublattice anisotropic Heisenberg 
antiferromagnet for spin 1/2 is then 


U = <H> =-NW) 
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m - tjF + 


=■ + m > 

J(0) 


E. coth 




(51) 


where Kj is defined by eqiiation (B5) or (B6) . The sum ^ Ej^ coth(j3Ej^/2) can be 
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evaluated with the help of the integrals of appendix B. Using equation (24) and convert- 
ing from a sum to integrals give 



where I 2 and II 2 are given by equations (B2) and (B4), respectively, and K 2 is de- 
fined by equation (B5) or (B6) . Substituting 



and 
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- F 


into equation (51) and simplifying resvilt in 
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Defining a reduced energy by u = U/NJ(0) gives 


(53) 
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Specifically, 
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Specifically, 


and 


u = + 7 ]F ( m + i ) (RPA) 


(55) 


u - - 


+ 77 F ( 21 x 1 ^ ^ (Callen) 


(56) 


The quantity u is plotted against the reduced temperature t in figures 3 and 4 for the 
RPA and Callen approximations, respectively, for t? = 0, 0.'5, and 1. 

For both decoupling approximations and 77 = 0, u = 0 at T^, This confirms the 
persistence of some order even above T^ due to the transverse coupling of spins . The 
intersections of the cirrves for the RPA appear to result primarily from the lesser weight 
given to the transverse coupling of spins. For the zero-field case it is not possible to 
determine from u alone which decoupling scheme gives a more accurate description of 
the internal energy. A calculation and comparison of the specific heats is more suitable 
for this purpose. 

Values of u at the Neel temperature are presented in table I. The gro\md-state 
energy u^ is the value of the internal energy at 0 K. These values are given in table n, 
tc^ether with the values for the sublattice magnetization at 0 K. The ground-state energy 
for the isotropic case (tj = 1) in the Callen approximation agrees exactly with the value of 
-0.293 obtained by Lee and Liu. They use a somewhat different and more complicated 

method to evaluate the longitudinal correlation f\mction ^^ 2 ^ which appears in <H) 

(ref. 3). 


NEEL TEMPERATURE 


The transition temperature for an antiferromagnet (Neel temperatime) is defined as 
that temperature at which the sublattice magnetization vanishes. With Q = Q'A and 
Q' = m/r, the reduced Neel temperatm^e Tj^ is determined from equations (27), (28), 
and (36) as 
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where A = 1 for the random phase and 1 - 2tjF for the Callen approximations . It is 
easily shown that 


(Q’A) coth (Q’A) = 1 

Q'A-0 

► 

(Q’A) csch (Q’A) = 1 
Q’A-0 J 


(58) 


Then 



lim J_ 
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The calculation for the Callen approximation also requires 
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Then equations (58) and (59) are used to determine 

The calculated values of for t? = 0, 0. 5, and 1 are presented in table in. For 
T 7 = 1, both the RPA and the CaUen Neel temperature values agree exactly with those ob- 
tained previously for the antiferromagnet (refs. 3 and 11). The value of = 0. 5 for 
77 = 0 is the well-known value for the Ising model. It is noted that the reduced Neel tem- 
peratTires for the antiferromagnet are identical to the reduced Curie temperatures for the 
ferromagnet for both the random -phase and Callen approximations. 


CONCLUDING REMARKS 


Our results for the sublattice magnetization and the Neel temperatures agree with 
results obtained previously from high- and low-temperature expansions. This suggests 
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that the method for evaliiating the lattice sums which was successfully applied to the 
study of the anisotropic Heisenberg ferromagnet can be applied with equal validity to the 
two-sublattice anisotropic Heisenberg antiferromagnet. The res\ilt is an analytic ex- 
pression for the sublattice magnetization that depends on the decoupling scheme. This 
expression is valid over the whole temperature range. 

For spin 1/2, an exact expression for the internal energy due to the ordering of 
spins (eq. (51)) has been derived for the two-sublattice antiferromagnet and is indepen- 
dent of any decouplii^ scheme or approximation. This expression was then evaluated 
for the bcc lattice using the Flax -Raich method, for both the RPA and Callen decouplii^ 
schemes. The result is an analytic expression for the internal energy valid over the 
whole temperatxire range. In the Callen approximation, the value of the groimd-state 
energy for the isotropic model (77 = 1) agrees exactly with the value obtained previously 
from a low-temperature expansion. The internal energy in the presence of an applied 
field can be obtained by essentially the same methods. The analytic expression for the 
internal energy can then be used to calculate the specific heat and many other thermody- 
namic quantities and hence is a starting point for formulating the thermodynamics of the 
two -sublattice antiferromagnet. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 8, 1971, 

129-02. 
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total number of lattice sites 

spin operators 

temperature 

internal energy 
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structure factor 
nearest-neighbor vector 
anisotropy parameter 
reduced temperatxire 
correlation function 
Green function 


18 



APPENDIX B 


EVALUATION OF INTEGRALS FOR BODY-CENTERED CUBIC LAHICE 

The method for evaliiating the integrals which appear in the calculations of the cor- 
relation functions in the random-phase and Callen approximations are outlined in this 
appendix. The methods are similar to those applied to the ferromagnet (refs. 4 to 6). 
The integrals are 
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with ^ = K 2 /Kj^. For the bcc lattice, the structvire factor is 
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The denominator can be expanded as 
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where 
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Also, 
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The series corresponds to a complete elliptic integral of the first kind (ref. 4): 




2n 


n=0 


where 




and /Cj^ and are conjugates. Then 
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The integral II ^ can be written in terms of Ij^ as 
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and is thereby easily evaluated. The evaluation of Ig proceeds as follows: 
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Then 
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where 



Using the same expansion as in the evalviation of , 
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Substituting 1/F^ = S^/(R^ +S^), where S = /3mKj^/2ff, gives 
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® (R^ + S^)^ 
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(B24) 


The sums over R are evaluated by the method of Laplace transforms (see refs . 4 to 6 
and references cited therein) . Sufficiently accurate results are obtained by keeping only 
the first few terms in the sum over n. Specifically, the sums over R needed to calcu- 
late Ig and II2 are 
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Finally, 


I = - i-/l + (1 + if.) coth (77S) + ifM.) csch^ (ttS) 

2 77S\ 8/ V 16/ 16 


(B27) 


The integral Eg is evaluated in a similar way. For the RPA, nS = m/r and k = rj. For 
the Callen approximation, jrS = m(l - 2 t 7 F)/t and ^ = (77 - 2F)/(1 - 277F). 
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APPENDIX C 


EVALUATION 


OF THE TERM 2 £ (sf,S| ) 

tg 


For spin 1/2, the term 2 ^ 2g^®lf®2^ evaluated exactly. The oper- 
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ator S can be written as 
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equation of motion given by equation (2) for a,fi- 1: 


EA 


IZ , Im 



- 2 



((^ 


qZ q-t- 

®2g^U 





'qZ q+ . 

^lZ^2g’ 



(C3) 


2 

Substitutir^ equations (Cl) on the right side for S gives 
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(C4) 


, Im - , 2g^\l , Im ^^2g, Im^ ■*■ ^ , 2g 

g g 

^ ^'^2g2gU , Im “ 11 2g, Im^ 

where ^2g2gll , Im “ ^^®2g®2g® IZ ’ ^ Green function of higher order. 

To simplify the calculation, the spectral density p for a Green function is 
duced. It is defined in reference 7 by 


p(A,B-,E) = 


i lim 
e-0 




The corresponding correlation functions are calculated from equation (5). For 
eqviation (5) can be written as 


/ ' OO 2 

- 1^ p(A,B;E)dE 

00 


Equation (C4) can be written in terms of spectral densities. Setting 1 = m = f 


Ep(lf,lf-,E) +^p(2g,lf;E)] + 2 

g ’ g ’ 


2g 


X [p(2g2glf,lf;E) - 77p(lflf2g,lf;E) 


where 


p(lf,lf-,E) = 


i lim 
e-0 


^If, If^® ^ " "^If, If^® " 


p(2g2glf,lf-,E) = 


i lim 
e-0 


^2g2glf, If " '^2g2glf , If^® “ 


intro- 

(C5) 

t = t’ 

(C6) 
results in 

(C7) 

(C8) 

(C9) 
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Multiplying eqioation (C7) by - 1) integratir^ over E from -«> to +°°, and 
summing over f give 


If 

f 


(e^E _ i)"^Ep(lf,lf;E)dE 2g 

f,g 





+ rj 



2g)j 

f,g 




(CIO) 


On the right side, jfS = 0, because for S = 1/2, (S“)^ =0 for any state 

Fourier transforming p(lf, lf;E) to reciprocal space yields 


Y, pdf, If; E) = Y Pii(k;E) 
f 


it 


Irt' 


E + ie) - ^ (k; E - ie) 


11 ' 


(Cll) 


Rearranging equation (CIO) gives 


‘Y - D'Vu^-.EjdE 

f’g k f,g 


2g 


X 


fefSlf) ^ (C'2) 


Then from equations (C2) and (C12), 


f,g f,g k 


xEp^^(k;E)dE (C13) 
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Also, 


and 

f,g 

The integral over E is evaluated by using the Green function of eqmtion (17) 
spectral density p^j^(k;E). The result is 

I. ~ \ ^ 

k k 


Substituting equations (C14), (Cl 5), and (C16) into (C13) gives 


f,g 


NJ(0) 

2 


m + 77 F + 




(C14) 


(C15) 


in the 


(C16) 


(C17) 
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TABLE I. - REDUCED INTERNAL ENERGY EVALUATED 
AT NEEL TEMPERATURE FOR SPIN 1/2 ANTI- 
FERROMAGNET WITH BODY -CENTERED 
CUBIC LATTICE 


Anisotropy 

Approximation 

factor, 

n 

Random phase 

Callen 


Reduced internal energy, u = (H)/NJ(0) 

0 

0 

0 

.5 

-.0043 

-.0097 

1.0 

-.0706 

-.0706 


TABLE n. - REDUCED INTERNAL ENERGY AND SUBLATTICE MAGNETIZATION 
EVALUATED AT T = 0 K FOR SPIN 1/2 ANTIFERROMAGNET 
WITH BODY -CENTERED CUBIC LATTICE 


Anisotropy 

factor, 

n 

Approximation 

Random phase 

Callen 

Sublattice 
magnetization , 
mo 

Reduced internal 
energy, 

Uq = <H>/NJ(0) 

Sublattice 

magnetization, 

mo 

Reduced internal 
energy , 

Uq = <H)/NJ(0) 

0 

0.5000 

-0. 2500 

0.5000 

-0.2500 

.5 

.4293 

-.2584 

.4907 

-.2597 

1.0 

.4706 

-.2906 

.4706 

- . 2934 
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TABLE in. - REDUCED NEEL TEMPERATURE 


FOR SPIN 1/2, TWO -SUBLATTICE, ANISO- 
TROPIC HEISENBERG ANTIFERROMAGNET 
WITH BODY-CENTERED CUBIC LATTICE 


Anisotropy 
factor , 

Approximation 

Random phase 

Callen 

Reduced Neel temperature 

0 

0.5000 

0.5000 

.5 

.4830 

.4970 

1 

.3589 

.4601 



Figure 1. - Magnetization as function of reduced 
temperature for random-phase approximation. 





.1 .2 .3 .4 .5 

Reduced temperature, t= kT/J(0) 


Figure 2. - Magnetization as function of reduced temperature for 
Callen approximation. 



0 .1 .2 .3 .4 .5 

Reduced temperature, t= RT/J(0) 


Figure 3. - Internal energy as function of reduced 
temperature for random -phase approximation. 




.1 .2 .3 .4 .5 

Reduced temperature, t= kT/J(0) 


Figure 4. -Internal energy as function of reduced 
temperature for Callen approximation. 
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